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Abstract : 



O 
^ . 

bJO, The onset of small-scale fluctuations around a steady convection pattern in a rotating baroclinic annulus filled with 

C/3 ' air is investigated using Direct Numerical Simulation. In previous laboratory experiments of baroclinic waves, 

_ ^ , such fiuctuations have been associated with afiow regime termed Structural Vacillation which is regarded as the 

C/3 ■ first step in the transition to fully-developed geostrophic turbulence. 



Resume : 



Le developpement de fiuctuations de petite echelle dans un ecoulement convectif permanent d'air en cavite tour- 
^ ' nante barocline est analyse par simulation numerique directe. Des travaux experimentaux anterieurs ont associe 

■"sj" \ de telle s fiuctuations au regime de vacillation structurelle, considere comme la premiere etape vers la turbulence 

C^^ ' geostrophique. 
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-^ , 1 Introduction 

X 

c^ ■ The transition to disordered behaviour in the form of 'Baroclinic Chaos' provides an impor- 

tant prototypical form of chaotic transition in fluid dynamics. This is of particular geophys- 
ical relevance in the context of understanding the origins of chaotic behaviour and limited 
predictability in the large-scale atmospher es of the Earth and ot h er terrestrial planets, such 
as Ma rs, and in the oceans (e.g Pierrehumb ert & Swanson I (1 1995b : iRead et all (11998); ReadJ 



(|2001|) ). For many years, the differentially-heated, rotating cylindrical annulus has proved a 



fruitful means of studying fully-developed, nonlinear baroclinic instability in the laboratory 



(|Fowlis & Hide I (|1965|) . lFruh & Read I (| 1997 ')). Transitions within the regular wave regime fol 



low canonical bifurcations to low-dimensional chaos, but disordered flow appears to emerge via 
a different mechanism involving small-scale secondary instabilities. Not only is the transition 
to geostrophic turbulence less well understood than those within regular waves, but also the 
classification and terminology for the weakly turbulent flows is rather vague. Various terms 
applied include Structural Vacillation or Tilted-Trough Vacillation which both refer to fluctu- 
ations by which the wave pattern appears to change its orientation or structure in a roughly 
periodic fashion. These fluctuations have been explained by the growth of higher order radial 
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Figure 1: Summary of the regimes obtained by the numerical model for a temperature difference of 
AT = 30K. The axis at the bottom shows the Taylor number while the axis at the top shows the ratio of 
the centrifugal term to gravity at the inner radius. 



mode baroclinic waves (Weng & Ba rcilon I (|1987h ). by barotropic instabilities, or by small-scale 
second ary instabilitie s or eddies, which lead to erratic modulations of the large-scale pattern 
(IHigne tt et a/.l (|1985h ). Subsequent development within this so-called 'transition zone' leads to 
the gradual and progressive breakdown of the initially regular wave patter n into a n incr easingly 
disordered flow, ultimately leading to a form of 'geostrophic turbulence' ([Read I (|200l|)). 



2 The numerical model 

The physical model comprises a body of air contained between two vertical, coaxial cylinders 
held at constant temperatures and two horizontal insulating rigid endplates separated by a dis- 
tance of d. The inner cylinder at radius r = a is cooled (T^) and the outer cylinder at radius 
r = 6 is heated (AT = T^ — Ta = 30 K). The whole cavity rotates at a uniform rate around the 
axis of the cylinders. The geometry is defined by an aspect ratio, A = d/{b — a) = 3.94 and a 
curvature pararneter, R r = {h + a)/{h — a) = 3.7, corresponding to the configuration used by 
Fowlis & Hid e ' ('1965') in their experiments with liquids. Following a long-established tradition 
(e.g. Fowlis & H ide (1965)), the analysis of the flow in a fixed geometry for a given fluid, here 
air at ambient temperature with Pr = 0.707, is made by varying two main parameters: the 
Taylor number Ta and the thermal Rossby or 'Hide' number 6, 



Ta 



4n\b 



p'^d 



e 



gdaAT 



The Navier-Stokes equations coupled with the energy equation via the Boussinesq approxima- 
tion are solved using a pseudo-spectral coUocation-Chebyshev Fourier method associated with 
a second order time scheme. 



3 The onset of SV regime 

Figure \\\ shows the equilibrated amplitude of the dominant azimuthal mode against the Taylor 
number. For steady solutions, the equilibrated amplitude of that wave mode is shown while 
for time-dependent solutions, the maximum and minimum as well as the mean amplitude are 
given. As such, the figure is a bifurcation diagram representing the sequence of three basic 
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Figure 2: Azimuth-height maps of eddy fields (with azimuthal flow removed) of a wave-3 simulation at 
(i) Ta = 0.235 x 10^ ; {ii) Ta = 4. x 10^ . (a) Pressure (contour interval = 0.25 dimensionless units), 

(b) Temperature (contour interval = 0.05, normalised to an imposed temperature difference of ±1.0) and 

(c) Vertical velocity (contour interval = 0.5 dimensionless units). 



solutions, the axi symmetric solution on the x-axis and the m = 2 and m = 3 solutions pro- 
jected onto the same plane. The instability of th e axisymmetric solution and the m = 2 solution 
branch are reproduced from the earlier study in lRandriamampianina et al\\200w . The present 
study focusses on evolution of the m = 3 solution branch when increasing progressively the 
rotation rate up to Ta = 5. x 10^. A transition in the flow structure is observed, where the 
three-dimensional structure of the flow responds to a shift in the balance between gravity and 
the centrifugal force, quantified by the local Froude number, Fr = Vt^r/g, and the dominant 
physical process changes from baroclinic instability to convection due to radial buoyancy. The 
transition of this convection to chaotic behaviour is fundamentally different from that observed 
in the transition to the chaotic flow observed at lower rotation rates. Rather than via a sequence 
of low-dimensional, quasi-periodic states, the large-scale convection developed small-scale in- 
stabilities, which have been previously suggested as the origin of Structural Vacillation (SV) 
within the transition to geostrophic turbulence. 

The Froude number taken at the inner radius Fva = Vl?a/ g was used as the secondary 
axis on top of the graph in Figure [U and the location of Fr^ = 1 is highlighted by the dotted 
vertical line at Ta = 414, 000. The first vertical line, at Ta = 239, 000, indicates where the 
centrifugal force equals gravity at the outer cylinder, Fr^, = 1. Below that Froude number or 
Taylor number, gravity is the stronger of the two forces everywhere in the annulus. It can be 
seen that Fr^ = 1 corresponds to a place on the solution branch where the radial structure of 
the steady wave changes. At the bifurcation, a higher radial mode of the same azimuthal wave 
grows to a finite amplitude while the first radial mode appears saturated around Ta ~ 10^ and 
even decays. In terms of the physical quantities, this development results in strong velocity and 
temperature gradients near the boundaries and the concentration of the radial transport of heat 
and momentum in narrow plumes or jets. At the highest values of Taylor numbers explored, 2 
to 5 X 10^, that flow then developed temporal fluctuations. 

This change from a regular steady wave dominated by baroclinic instability to the transition 
zone characterized by small-scale fluctuations can be seen from the variations of eddy variables 
at two different values of rotation rates in the azimuth-height maps at mid-radius in Figure [21 
At Ta = 0.235 x 10^, the solution shows a pressure field which tilts westward with height and 
the temperature and vertical velocity fields which tilt eastward in such a way that the strongest 
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Figure 3: Maps of the odd EOFs of the first six pairs EOFs of temperature at mid-height from a simulation 
of baroclinic waves in air at Ta = 4.0 x 10^. EOFs are shown in order of decreasing contribution to the 
total temperature variance at this level. 



upwelling is located around the strongest temperature gradient and strongest pressure gradient 
near the lower boundary. This behaviour is that of a typical baroclinic steady wave. At higher 
rotation rate Ta = 4. x 10^, the pressure field clearly shows very little phase tilt with height 
while the temperature field exhibits highly concentrated plumes of hot and cold air which have 
very little slope with height, interspersed with broad regions with relatively very weak hor- 
izontal thermal gradients. The thermal plumes are aligned azimuthally with the regions of 
strong azimuthal pressure gradient, consistent with optimising the correlation between radial 
(geostrophic) velocity and temperature perturbations. Similarly, the regions of strong upward 
vertical velocity are also concentrated close to the strongest positive temperature anomalies, 
and vice versa, consistent with an optimisation of w'T' (where w' and T' represent departures 
in vertical velocity and temperature from their azimuthal mean values). Eventually, the decrease 
of the tilt results in a virtually vertical, 'barotropic' structure. Regions of strong downward mo- 
tion (necessary to satisfy mass conservation) are concentrated in plumes or jets adjacent to the 
strong upwelling jets, forming 'cross-frontal' circulations which show some similariti es with 
those inferred for atmospheric frontal regions in developing cyclones (|Hoskins I (|1982h ). The 
overall impression is that the flow in this region of parameter space is strongly nonlinear and 
much modified from the simple, linearly unstable Eady solution found at much lower Taylor 
number. 



4 The SV characteristics 



The spatial structure associated especially with the onset of structural vacillation is of particular 
interest, since it may be associated with bifurcations involving quite different spatial modes. As 
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Figure 4: Power spectra from the time series of the temperature after removing the steady wave: (a) near 
the inner wall, (b) at mid-radius, and (c) near the outer wall. All show a line ~ f~^ but (c) also shows a 
line ~ /~^' ^ at the lower frequencies. 



a means of isolating the dominant pat terns in the flow, we h ave made use of a form of Empirical 
Orthogonal Function (EOF) analysis (jPreisendorfer I (|l988|) '). in which time sequences of spatial 
maps e.g. of temperature were analysed to obtain the covariance matrix of every spatial point 
with every other in the two-dimensional field. 

Figures |3] show the first six pairs of EOFs at Ta = 4.0 x 10^, which together account for 
a little over 98% of the temperature variance. However, the first six EOFs appear as three 
conjugate pairs and represent 93.2% of the variance, indicating a slightly broader spread into 
the higher order patterns. The first two EOFs, e.g. Fig. [3]; a), correspond to pattern dominated 
by m = 3, where the radial structure seems somewhat more concentrated towards the side 
boundaries, with only weak amplitude throughout most of the interior radii. EOFs 3 and 4, 
Fig. Sb), show a pattern dominated by m = 6, but with even more contorted phase shifts with 
radius. The overall impression is of a pattern represented by m = 6 in azimuth and / = 4 in 
radius (where / is a nondimensional radial wavenumber). Subsequent EOFs display patterns 
with a reasonably simple / = 1 radial structure (except for (c), which also shows evidence for 
I = 4 dependence on r in EOFs 3 and 4), and increasing azimuthal harmonics of m = 3. 
The higher order EOFs in Figs [3tc) and (d) generally have a much simpler radial structure, 
with a simple maximum in amplitude near mid-radius and only weak phase-tilts with radius. 
They again appear in conjugate pairs, and represent other azimuthal harmonics of m = 3. The 
complex phase variations with radius in EOFs 3 and 4 demonstrate the growth of higher radial 
modes as previously supposed in Figure [H though the patterns in (c) and (d) of Figs [3] would 
seem to suggest such anti-correlations ought to be strongest for m = 6 rather than the dominant 
m = 3. 

Spatially averaged power spectra from the temperature residuals, after removing the steady 
wave, at mid-height and three radial positions, at mid-radius and around 15% from each wall, 
are shown in Figure IH The main frequency (~ 0.03T^\ with the time scale T = l/2i7) is 
that of the emission of the perturbations described by EOFs 3 and 4. The spectrum falls off 
at all radial positions, with some distinctive peaks over the general decay. While the decay at 
higher frequencies is largely consistent with a f~^ law at all postions, a slight flattening in the 
frequency range between that of the main perturbation and about O.IT^^ can be observed near 
the outer wall. The behaviour of the spectrum in that range appears to be closer to a /^^^'^ law. 

The spectral evidence suggests that the flow investigated here consists of a fairly steady 
large-scale convection pattern in the form of three pairs of hot and cold radial jets. From those 
jets, smaller perturbations are emitted at relatively regular intervals where a jet approaches a 
wall. The overall flow field responds in a cascade of faster (and presumably smaller) fluctuations 
which appear consistent with two-dimensional, quasi-geostrophic turbulence over a wide range 
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of frequencies similar to results of DNS of geostrophic turbule nce. Recently, similar energy 



spectra were observed fo r geostrophic turbulence in a square box (|Lindborg & Alvelius I (|2000l) 



Waite & Bartello I (l2006h ). 



5 Conclusions 

In the present work, we find the emergence of a new steady wave solution with complex radial 
structure, associated with a change in the balance of forces to a centrifugally-dominated regime, 
followed (at high Taylor number) by a further transition to a form of structural vacillation asso- 
ciated with regular m = 3 eddies driven by centrifugal buoyancy. The computations have been 
carried out on the NEC SX-5 of the IDRIS (CNRS, Orsay, France). The authors are grateful 
to the British Council and the CNRS for funding the collaboration between the French and UK 
partners in a joint programme Alliance. 
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